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CP asymmetry, branching ratios, and isospin breaking effects of B ^ K*'y with the 

perturbative QCD approach 
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The main contribution to the radiative B — > 7^*7 mode is from penguin operators which are 
quantum corrections. Thus, this mode may be useful in the search for physics beyond the standard 
modeL In this paper, we compute the branching ratio, direct CP asymmetry, and isospin breaking 
effects within the standard model in the framework of perturbative QCD, and discuss how new 
physics might show up in this decay. 
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I. INTRODUCTION 

The large CP violation in B -^ J jijjKs decay mode 
predicted by the standard model with Kobayashi- 
Masukawa (KM) scheme has been verified by B factories 
at KEK ( High Energy Accelerator Research Organiza- 
tion) and Stanford Linear Accelerator Center (SLAC). 
The standard model predicts the CP asymmetries for 
B -^ J/ipKs and B -^ ipKg to be equal to sin 20i . How- 
ever, recent experimental data from Belle showed that 
these asymmetries differ by nearly 2ct; the averaged 
sin 2(j>i from Bell and BaBar in B ^ J/ipKg system 
is sin 201 = 0.736 ±0.049 [J and in B ^ ^isT, decay 
mode is sin2(/)i = 0.06 ± 0.33 ± 0.09 from Belle 0, and 
sin 201 = 0.50± 0.25to;ol from BaBar 3]. Experimental 
error is still large, so the situation is inconclusive, but if 
this result continues to hold, it implies existence of new 
physics beyond the standard model. 

In this paper, we want to concentrate on _B — > K*j 
decay mode. The decay mode has large a branching ratio, 
so the experimental error on the CP asymmetry has been 
getting small and is now down to several percent. 



Br{B" -^ K*°-f) 



Br{B^ -^ K*^j) 



iCP 



f (4.01 ± 0.21 ± 0.17) X 10-^ [4j 
|(3.92 ± 0.20 ± 0.24) X 10-5 [5] 



' (4.25 ± 0.31 ± 0.24) X 10-5 [4] 
(3.87 ± 0.28 ± 0.26) X 10-5 [5j 



-0.015 ±0.044 ±0.012 [4] 
-0.013 ±0.036 ±0.010 [5] 
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CP Asymmetry is defined as 



A -^^^ 



K*j) ~ V[B ■ 



^K*j) + r{B 






(1) 



and in general, theoretical predictions of the CP asym- 
metries depend less on hadronic parameters than those of 
branching ratios as many uncertainties cancel in the ra- 
tio. So comparing predictions for CP asymmetries within 
the standard model with experimental data may be ef- 
fective way to search for new physics. Many authors 
have pointed out for some time, that the CP asymmetry 
in this mode is very small. We can easily understand 
why the asymmetry is so small. In order to generate 
CP asymmetry, at least two amplitudes with nonvan- 
ishing relative weak and strong phases must interfere. 
This decay is mainly caused by an Oj-y operator, and 
other contributions which interfere with this contribution 
are small and the Cabibbo-Kobayashi-Masukawa quark- 
mixing matrix (CKM) unitary triangle is crushed, mak- 
ing the CP asymmetry in this decay mode very small. 
However, if we were to look for new physics, we need to 
be able to give a quantitative estimate of the standard 
model contribution to the CP asymmetry. For this pur- 
pose, we include small contributions which interfere with 
O7-),, including also the long distance contributions, for 
example, B -^ K*J/'ip -^ K*j. 

Furthermore, the isospin breaking effect Ao+ is also 
very interesting because it's size and sign are sensitive to 
the existence of physics beyond the standard model. 



An+ = 



r{B° -^ X*%) - T{B+ -> K*+j) 
r(B" -^ K*^-i) + T{B+ -^ K*+j) 
{TB+/TBo)Br{B° -^ K*°j) - Br{B^ 



K*+-f) 



(TB+/TBo)Br{BO -^ K*Oj) + Br{B^ 



(2) 



In order to test the standard model, we need to 
know if the penguin contribution within the standard 
model can explain the experimental data. Experi- 
ments show Ao+ = +0.012 ± 0.044 ± 0.026 in Bell ^ and 



Ao- = 0.050 ± 0.045 ± 0.028 ± 0.024 in BaBar [|. More 
precise data will become available in the near future, 
so the theoretical prediction of it's size and sign of this 
asymmetry should be pinned down. 

In this paper, we calculate the branching ratio, direct 
CP asymmetry, and isospin breaking effects in _B ^ K*^ 
decay mode, based on the standard model. First, we 
briefly review the concept of the pQCD in Sect^I and in 
Sect lllll we show the effective Hamiltonian which causes 
B -^ K*^ decay. Then we present the factorization for- 
mulas for the B -^ K*j decay mode in Sect II VI and in 
Sect0 we mention about the long distance contribu- 
tions. Next we will show the numerical results in Sect I VII 
and Sect lVIll is our conclusion. Finally in Appendix A, 
we present a brief review of pQCD. 



II. OUTLINE OF PQCD 

Theoretically, it is easy to analyze the inclusive B me- 
son decay like B -^ Xs^ because we can estimate the de- 
cay width, for example, by inserting the complete set for 
all possible intermediate states. The experimental and 
theoretical branching ratio oi B ^ Xg'y are 



Br{B - 
Br(B 






(3.52l°-jg) X 10-^ 
(3.57 ± 0.30) X 10"^ 



and this good agreement strongly constraints new physics 
parameters. However, inclusive decays are experimen- 
tally difficult to analyze because all B -^ Xg^ candidates 
should be counted. If we can directly calculate the ex- 
clusive decay mode B — j K*^, we ought to obtain many 
interesting results to test the standard model or to search 
for new physics. 

Perturbative QCD is one of the theoretical instrument 
for handling the exclusive decay modes. The concept 
of pQCD is the factorization between soft and hard dy- 
namics. In order to physically understand the pQCD 
approach, we consider B^ meson decays into K*'^ me- 
son and 7 in the rest frame of the B^ meson (Fig^ . The 
heavy b quark which has most of B^ meson mass is nearly 
static in this frame and the other quark, which forms the 
B^ meson together with the h quark, called the spectator 
quark, carries momentum of order 0(A) = 0{Mb — rrib). 
This b quark decays into the light s quark and 7 through 
the electromagnetic penguin operator and the decay 
products dash away back-to-back, with momentum of 
0{Mb/2). (This process is depicted in Figda).) K*° 
meson is composed of s quark and a spectator quark. In 
order for the fast moving s quark and slow moving spec- 
tator d quark to form a K*'^ meson and nothing else, the 
spectator quark must be kicked by the gluon, so that the 
s and d quark have more or less parallel momenta in the 
direction of K*^. (This process is depicted in Fig^Jb).) 
Since the invariant-mass square of this gluon is the order 
of 0{AMb), we can treat this decay process perturba- 
tively. 




s^^K 




■^*o 



FIG. 1: The left figure is no gluon exchange diagram, s and 
spectator d are not lines up to form an energetic K* meson. In 
order to hadronize K* meson, one gluon with large q^ should 
be exchanged. 




FIG. 2: Example of annihilation diagram which produces 
strong phase through the branch-cut. 



There is also diagram shown in Fig [3 This can 
also be computed in the pQCD approach. The dia- 
gram can be cut along the dotted line indicating the 
presence of the physical intermediate state. This re- 
sults in a strong interaction phase which can be com- 
puted. The direct CP asymmetry is caused by interfer- 
ing some amplitudes which have relative weak and strong 
phase, and it can be written in the form proportioning 
to sin {0wi — Ow2) sin {5gi — 5s2)'- in short, it depends on 
both weak and strong phases. 
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We can determine the strong phases by using the 
pQCD approach, then we can extract the information 
about the weak phases and examine the standard model. 
A more detailed review for the pQCD approach is in Ap- 
pendix ^ 



III. KINEMATICS FOR B -^ K*'y DECAY MODE 

The effective Hamiltonian which induces flavor- 
changing b ^ sj transition is given by (g| 



H, 



eff 



^-^ E v,,v;g[c^{t,)o['\t,) + C2ifi)o^^\t.) 



V2L 



q—u,c 



h.c. , 



(3) 



Oi = {siqj)v-A{qjbi)v~A, 

02 = {siqi)v~A{q]bj)v^A, 

03 = {sibi)v-A^{qjqj) 



V-A, 



04 = {sibj)v-A^iq]qi) 

q 

05 = {sibi)v-A^{qjqj) 

q 

Oq = {sibj)v-A^{q]qi) 



V-A, 



V+A, 



(4) 



V+A, 



then the b and s quark momenta are pb = Pi — fci and 
Ps = P2 — k2, and we neglect the masses of the Ught 
quarks and identify the b quark mass with the B meson 
mass in calculations of the hard scattering amplitudes. 
The term proportional to Af, = AIb — rub is generated by 
higher order effects, so we included this effect in our error 
estimate. 

From here, we extract the formulas for decay ampli- 
tudes caused by each operators, 

M = {F\ H^ff I /) 

= ^E^C7i.MC.(M)(i^|0.(/i)|/) (11) 



O77 = —^Siu'^'' {msPL +mbPR)biFf_i^, 
Osg = ■£^-S:<j'^''{m,PL+mbPR)Ttjb,Gl,, 

where P^ = (1 =F 7^)/2. We define the B meson and the 
K* meson momenta Pi and P2 in the light-cone coordi- 
nates 

P= (p^,P",Pt) = ( — -=— , — -j^,{p^,p^)\ (5) 

within the B meson rest frame as 

Fi = (P+,Pf,PiT) = ^(1,1, Ot), (6) 

F2 = (P+,^2^^2t) = ^(0,1,0t), (7) 

and photon and the K* meson transverse polarization 
vectors as 



1 



e;(±) = (0,0,^(Tl,-*) 
e*i,.{±) = (0,0,i=(±l,-z) 



(8) 



and they can be decomposed into scalar and pseudoscalar 
components as 

M = (e; • e*K,)M^ + ^e^,+^e;^eJ^^.^f^. (12) 



IV. FORMULAS 

In this section, we want to show the explicit formu- 
las of the decay amplitudes caused by operators given in 
SectHnl 



A. Oj-y contribution 

If we define the common factor as 

i.(0) ^ ^lv:,yr^Ml 



%/2 7r 



V^b^cs'-^F^y^B^ 



(13) 



where Cp is color factor, and S,i as V*^y^s/V*fycs^ the 
decay amplitude A/7-y in Fig|31can be expressed as follows. 



Throughout this paper, we keep only terms of or- 
der rK* in the computation of the numerator, where 

TK' =Mk'/Mb. 

The fractions of the momenta which have the spec- 
tator quarks in B and K* meson are Xi = kf /P^ and 
X2 = fc^ / P2 I so the momenta of these spectator quarks 
are expressed as follows, 

, _ ^ Mb 

ki = {k{,k^ ,kiT) = {—^xi,0,kiT), (9) 

^2 = ik:^,k2,k2T) = iO,—7=-X2,k2T), (10) 



(a) 



(b) 



Oj, 



FIG. 3: Feynman diagrams of electromagnetic penguin oper- 
ator Or-/ . A photon is emitted through the Or-, operator, and 
one hard gluon exchange is needed to form hadrons. 



m/^("' = -M^^^"^ 



= -2F(")Ct / dxidx2 / 6id6i52d62</'i3(2:i,6i)5*(2;i)a,(i?)e[-^«(*?)-'5^*(*?)lC7(i?) 



M. 



Sib) 

77 



X i/("' (A762, Brbi,Brb2) tk- Wk- {^2) + <i>K' {^2)] (i? = max(A7, B7, 1/61, 1/62)) 
= -2F^"^^t I dxidx2 fhdh b2db2<f>Bixi,bi) Stix2) a,(4)e[-^«(*')~^^*(*')]C7(4) 

X H^'^ (A7&1, C761, C7&2) [(1 ~ 2x2)rK^ {Cbl, {X2) + rK^ (X2)) + (1 + X2)^l, {X2) 

(4 = max(A7,C7,l/6i,l/62)) 



ij(") (A762, 5761,5762) = A'o(A762) e(6i " 62)A'o (5761) /o (5762) + d{b2 - bi)Kn (B762) 4 (5761) 



r(fc) 



(a)/ 



H)">{Ajbi,C7bi, C762) = 5) ^^761, C761, C762) 
A? = X1X2MI, B^ = a;iA/B' ^77^ ^ x2Af| 



Here Kq, Iq are modified Bessel functions wliich come 
from propagator integrations. The meson wave func- 
tions are not calculable because of its nonperturbative 
feature. But they should be universal since they absorb 
long-distance dynamics, so we can use the meson wave 
functions determined by some approaches. We use in this 
paper a model B meson wave function which is shown to 
give adequate form factors for B -^ Kn decays ^ ilOi] . 
and K* meson determined by light-cone QCD sum rule 
[ill Il2| . Their explicit formulas are shown in Appendix 
B. 



J 



(14) 



(15) 
(16) 

(17) 
(18) 



B. Osg contribution 



Similarly, we can calculate the Ogg contributions as 
follows. In these cases, a hard gluon is emitted through 
the Osg operator and glued to the spectator quark 
line (Fig0J. In the following formulas, Qq express the 
electric charge of the external quark: Qu — 2/3 and 

Qd = Qs = Qb = -1/3. 



M 



S{a) 



-M 



P(a) 



= -F^^ktQb I dx^dx2 fhdb, b2db2 <^fl(si,foi)5',(xi)Q,(ig)e["^«(*8)-^^*(*«']C8(i|)i:f^'"'(yl8fe2,B8&l,B862) 

X [X1(J>K*{X2) +rK*X2{(l)K'{x2) + (t>'K'{^2))] (tg = max(A8, -Bs , 1/bi , I/&2)) (19) 

= -F'^'CiQ. / dxidx2 I bidbi b2db2 (Peixi, bi)St{x2)as{tl) e[-^«<*«'-®^*<*«'lc8(t8)^8''(^86i, Cg&i, C862) 

X [-3x2rK'{(l>K-'(x2) + (l)K'{x2)) + {2x2-xi)cj,Jc*(x2)] (tg = max(A8, C8, I/61 , 1/62)) (20) 



Mf<'=' = -A/> 



= -F 



P(c) 



itQ, I dXidx2 fbidbi b2db2 <t>B{xi,bi)Stixi)asitl) e[-^«('B'-^^*'*«%8(i8)^8''(\/l^62,£'8&l,F'8&2) 

X l-xicj^Jc-. {X2) + X2rK' (<t>K' {X2) + 4>K' {x2))\ Ul = max(y^l A'i 1, Ds, 1/61, llb2)\ (21) 



(a) 



(b) 



(c) 

5 - 



5 - 



FIG. 4: Feynman diagrams of chromomagnetic penguin operator Oss- A hard gluon is emitted through the Osg operator and 
glued to the spectator quark hne. Then photon is emitted by bremsstrahlung of external quark lines. 



M^W ^ ^F^°^^,Q, 



I dXjdX2 fbidbi b2db2 cl,Bixi,bl)St{x2)as{tt) et'^^ ^'^ )-Sk- (tg)] Cg(id) J/W (^j^^^ ^ ^^^^^ ^^^^ 



X [ 6x2rK* (pK' {^2) + (2 + a;2 - xi)(j>'^, (X2)] 



j^P(d) ^ F^^'^tQq I dxidx2 [ bidbi bidbi 4>Bixi,bi)Stix2)as{tt) e[-^«''8'-^'^*'*«']c8(i8)i?8'''(\/i^&i. -Bsbi, 



^862 



X [ (2 + X2 -Xi 



. (X2) + Gx2rK'<f}K' (2:2)] (4 = miix{^\A'g^\,Es, 1/fei, 1/62)] 



Hi''\Asb2,Bsbi,Bsb2) = 7^0(^862) [0(61 - b2)Ko{Bsbi)IoiBsb2) + (61 ^ 62 
Hi''\Asbi,Csbi,Csb2) = YKo{A»bi)\^9{bi-b2)H^o'\Csbi)MC8b2) + (61 ^ 62) 
Hi"\^/\Ai\b2,D8bi,Dsb2) = 9{A'i) Ko{^/\Ai\b2)[e{b^ - b2)Ko{D8b^)h{Dsb2) + (61 ^ 62) 

+e{-A'i) i'^Hl^\.J\A^\b2)[e{b^ - b2)Ko{D8bi)Io{Dsb2) + (61 ^ 62 



Hi''\,/\Ai\b^,Esb,,Esb2) = e{A'i) i^Ko{^/\Ai\b^) [e(6i - b2)H^o'\Esb^)Jo{E8b2) + (fei ^ 62) 

-ei-A'i) (|)'i/,f^(y^l6i)[e(6i-&2)^^'>(i586i)Jo(i?8&2) + (&i^&2) 

Al = xiX2M%, b1^mI{1 + xi), cI=mI{1-X2), A'i ^ {xi ~ x2)mI, dI ^ xiMb^ , EJ ^ x2Mb^ 



(22) 

2) 
(23) 

(24) 
(25) 

(26) 

(27) 
(28) 



C. Loop contributions 

1. Quark line photon emission 

Next we want to mention about charm and up quark penguin contributions (FiglFJ. The subtitle like "Quark line photon 
emission" means that a photon is emitted through the external quark lines. We define the c and u loop function in order that 
the b —> sg vertex can be expressed as s^'^{l — 'y'')Ifj.vb. It has the gauge invariant form 13] and the explicit formula is as 
follows, 



^ ^Liy — ~^ 2"^ ^M^^ f^^f^u) j uXXyL X) 



Jo 



9T\,2 



87r2 



I."' 5'm'-' f^fif^u) 



G{m,,k ,t) 



1+log 
2" 



TTlf — x{l — x)k^ 



i2 



(29) 







FIG. 5: Feynman diagrams of "Quark line photon emission". The charm or up loop go to gluon and attach to the spectator 
quark hne. A photon is emitted through the external quark lines. 



where k is the gluon momentum and irn is the loop internal quark mass 

r^i 2 ,2 ,x n, ,2N2r5 , 4m^ 

G(mi ,k,t) — 0{—k ) 



4:m1 , V'l-4m2/fc2 - 1 



2 1-5 4m- ml 2mi\ 4m- ^l 

+e(fc')e(4m?-fc')|[ 

^a(iJ A 2^2^5 _L 4m? m? / 



2r5 Ami , mf „A 2mt 



4m| 
fc2 



1 arctan 



\/4m2/fc2 - 1 



! + ■ 



4m| 
fc2 



In ^ . '/ — + JTT 



1 + ^1 - 4m2/fc2 



(30) 



The loop function G has the dependence of gluon momentum square of k^ . But there is no singularity when we take the limit 
of fc -^ 0, so we can neglect kr components of k^ in the loop function G. 

Then the "Quark line photon emission" contributions can be expressed as follows. 

X G(mi, --^2,42) - -\xiX2rK'{4>K*{x2) ~ (f>'K'{x2)) (^2 = max(A2,-B2, 1/&1, I/&2)) 



(31) 



Mfl"^ 



^M{ 



p(b) 



_Q±p{0) ^^ j rf2;^da;2 Kidbi62d&20s(a;i,bi)C'2(4)a,(4)St(a;2)e[-^«('^'-^^*('^']Hf'(A2bi, 6*261, C2b2) 
><\G{mi,-A%,t\) - -UxirK*{(f>'ii'{x2) + (t>'K'{x2)) + 'ixiX2(f>]i'{x2)\ ft| = max(A2,C2, 1/61,1/62)) (32) 



Ml 



S(c) 



-M-(=) 



= ^f'°) e. / da;irfa;2 /6id6i62d62<^s(a:i,6i)C'2(t2)a.(t2)5't(a:i)e[~^«(*2)-SK.('2%W(y^^2^^^^^^^^^2^ 
x[G{mi,~A2\t'2) ~'^[x2rK'{fK'{x2) + <t>K'{x2)) - xi<t)l,{x2)\ (tl = ms.^{^f\Aj'\ , D2 , 1/61,1/62)) (33) 



(d) _ Qq p(0) 



M^^"" = 



^f(») ^i I da;id:c2 /" 61 d6i62d620s(:ri,6i)C2(t2)a.(t2)Si(a^2)e[-^«(^''-^^*(''%^"H^A^6i, £261,^262) 
X [G(m?, ~Ai, 4) - |] [xarK* [3(1 + Xi^K-- {X2) - (1 - X2)4>k' (2^2)] + 3(a;2 - a;i)0^. (0:2)] (34) 



(a) 



(b) 





FIG. 6: The Feynman diagrams of "Loop line photon emission". The charm or up loop go to gluon and attach to the spectator 
quark line. A photon is emitted by bremsstrahlung of loop quark line. 



X G{m'^, -A2 ,12) - n X2rK* [-(1 - X2)(t>"K* {X2) + 3(1 + X2)(I)K' {^2}] + S{X2 - Xi)(j>K* {X2) 



ta = max(J|A22|,S2, 1/61,1/62; 



(35) 



H'2'"'(A2b2,B2bi,B2b2) 



H^o''\Asb2,Bsbi,Bsb2) 



'-8 



Hl2''>{A2bi,C2bi,C2b2) = H^''>{Asbi,C8bi,C8b2) 
Hi''\^J\A^\b2,D2bx,D2b2) = H''^\\J\^\b2,DsbuDsb2) 
Hi''\,/^\bi,E2bi,E2b2) = Hi''\,/\A^\bi,Esbi,Esb2) 



(36) 



Al=xiX2Ml,, Bl = {l + xi)Ml, Cl = {l-X2)Ml, Ai = {x^ - X2)mI, dI = xiM%, eI = X2MI (37) 



2. Loop line photon emission 

Next we consider the "Loop line photon emission": a photon is emitted through the c or ii quark loop line. We sum up 
Fig|SJa) and|^b), the 6 -^ sg-y decay amplitude is expressed as 



A{b -> 537) = e'ii[q)ea{k) s{p')I^^b{p), 
where vertex function /^„ is defined as follows p^HsHlql. 

I'lliy = Fi {k ■ q)e^^p„ {q'' - k'')'y'' + q^eij,p„s q''k''^^ - k^e^p^g q''k'^j^ 
+F2 k„epp„s q'^k'^-y + k^ ep^pa q'^'y'^ L, 



(38) 



(39) 



Fi = -'^egT^Gp C dx [' ^ dy 



F^ = 



37r2 

i2y/2 
'1^ 



xy 



ml — 2xy{k ■ q) — k^x{l — x) ' 



egT"-GF dx dy 



x{l — x) 



mq — 2xy{k ■ q) — k'^x{l — x) 



(40) 
(41) 



where L = (1 — 7'')/2, k is the gluon momentum and q is the photon one. Then the amplitudes can be expressed as follows: 



Ml = -^^'"^^'Z '^^Z <^y I dxidx2 fbidbi(l>B{xi,bi)C2{t2)a.{t2)e^-^'''^''^-^'''^''^^H2{biA,bi^/\B^\ 

^ ir72 2 rJ'^2 (1 -2X2)TK'(t>KA^2) - (1 + 2xi)(t>K, {x2) + TK' (I>k* {X2)\ 

xyX2Mg — m^ L L J 

+x{l - x)\xlrK'{(l>K*{x2) + (I>K*{X2)) + 3xiX2(l>Jc* {X2) 



(42) 




FIG. 7: The energy scale of the loop contribution is t = niax(y4, ^/jSP, 1/bi), where A^ is gluon momentum , B^ is the loop 
momentum, and 61 is the transverse interval between the quark and antiquark pair in the B meson. 



M,^ = 



^F^°^^. [ dx f dy f dxidx2 /'bidbi<^i3(a;i,6i)C2(t2)Q.(i2)e'-''«(*=^'-^^-<'^"ii-2(6iA,foiV|B^) 
■J Jo Jo Jo J 

■\xyX2\rK-(l>K*{x2) - (1 + 2xi)(f)Jc*ix2) + (1 - 2x2)rK*(t)K'ix2)\ 

+x(l~ x) [xlrK* (<Pk* {x2) + (pK- (2^2)) + 3x1X20^* (3^2)] 
t2 = max(A, yi^. V^i)) 



xyx2M'^ — m? 



(43) 



^2(6171,6171521) = Ko{hiA)-Ko{b 



Ko(biA)^i-Ho{b 



B2|) (B^ > 0) 
B2|) (B^<0) 




(44) 



A' =XiX2Ml, 



B' 



■ X1X2MB 



-X2MI + 



^(1 



(45) 



In general, it is hard to estimate the u loop contributions accurately because of the nonperturbative hadronic uncertainties. 
For k^ < IGeV, nonperturbative correction to the u quark loop shown in Fig|7|is large and in fact uu pair may be better 
represented by resonances. On the other hand if k^ is large, the perturbative computation is expected to be reliable. 

In the pQCD approach, the factorization energy scale t is determined at each point of the integration, i.e. for each point 
{xi,X2,hi,b2). Then these variables are integrated over the entire physical region. So for each point (xi, 2:2, 61, 62), Qs(i) can 
be determined. Thus we can observe the contribution to the amplitude as a function of ctsit). Figures |S] and |U] show the 
distribution of aa{t) for a diagram with the c quark loop, and u quark loop, respectively. We can see that the major part of the 
c loop contribution comes from a perturbative region, on the other hand u loop contribution includes also a nonperturbative 
region. Since M2U and A/^ gets considerable contributions from the nonperturbative region, we introduce 100% theoretical 
error for these amplitudes. 



D. Annihilation contributions 



Tree annihilation 



We now discuss the annihilation contributions caused 
by Oi and O2 operators. The diagrams are shown in 
FigEl The operators Oi , O2 can be rewritten as 



01 = {SiUj)v-A{Ujbi)v-A = {Sibi)v-A{ujUj)v-A, (46) 

02 = {SiUi)v-A{ujbj)v-A = {Sibj)v-A{UjUi)v-A- (47) 



These annihilation contributions are tree processes: no 
hard gluons are needed because they are four Fermi inter- 
action processes and do not include spectator quarks which 
should be line up to form hadrons. However, these con- 
tributions are small because it has {V — A) ® {V — A) ver- 
tex: gets chiral suppression, and its' CKM factor is V^f,Vus, 
0{}?) suppression compared to V'^Vts and V*jyas- Defining 
a2{t) = C2{t) + Ci(t)/3, the each decay amplitudes are as fol- 
lows: 



jS{a) 



M^'-"'' = M2 



P(a) 



= -F'"'^u 



(())/ 



SV^QbfK'TT 



4:Ml 



tk* / dxi / bidb 



a2(t:)S't(a:i)e[-^«('^'Vs(2;i,6i)A'o(biA„) 
(i: = max(A„l/6i)) (48) 



M. 



S{b) 



F^'ku 



3\/6Qs/Bvr 



4M| 



-Vk* / dX2 j 



b2db2 



a2(tl)St{x2)e^-'^'<'^<'>\i'^H^^\b2Ba) 

X \{2- X2)lt>"K*{x2)+X2(l>°K'{x2) 



(49) 



M. 



p{b) 



TK* / dX2 I b2db2 



_ p{Q)f 3\/6Qs/B7r 
" ^" AMI 

a2(ta)St(x2)e[-^'^*('''ljjH«(&2Ba) 

X [x2(i)K' {3:2) + {2 — X2)(J)'k* {3:2)] 



ia = max(Ba, 1/62 



(50) 



FIG. 8: The horizontal hne is aa{t)/n and the vertical axis is the contribution from the each energy region in Os to the total 
decay amplitude. The left figure is the real part and right one is imaginary part of the c loop contribution. These figures show 
that we can compute the c quark loop contribution perturbatively. 



FIG. 9: The horizontal line is as{t)/Tr and the vertical axis is the distribution from the each energy region in Us to the total 
decay amplitude. The left figure is the real part and right one is imaginary part of the u loop contribution. These figures show 
that nonperturbative contributions are important, so we cannot accurately compute the u quark loop contribution and have to 
take into account the hadronic uncertainty. 




u u 




M. 



S{d) _ p(0)^ SVeQufBTT 



dxo 



62d&2 



X [{l+X2)rK'{x2)-{l-X2)rK'ix2)] (52) 



(o) 




(d) 




M. 



P(d) 



p(0)^ 'W&QufBjl 



dx2 I b2db2 



a2{ti)Stix2)el-'''^'^'i^h^H^o'\b2D,) 

X [{1-X2Wk'{x2)-{1+X2)(^''k*{x2)] 



ta = max(Da, 1/62 



(53) 



FIG. 10: Annihilation diagrams caused by Oi, O2 operators. 
In these cases, no hard gluons are needed because they are 
four Fermi interactions and do not include spectator quarks 
which should be line up to form hadrons. 



Al = {l+xr)Ml, 






--X2MI 



(1 - X2)Ml 



(54) 



2. QCD Penguin 



M^^"^ = -M^^"^ 



= F(°)e 



SV&QufK'TT 



4M| 



dxi / bidbi 



a2(t:)5't(a;i)e[-^«(*°']0i3(xi,bi)A'o(biCa) 
(i^ = max(Ca,l/6i)) (51) 



Next we mention the QCD penguin annihilation caused 
by O3 ~ Oe operators like in Fig llll Here we define 
a4(t) = C4(i) + C3(t)/3, aeit)^Ca{t) + Csit)/3. O3, O4 
have the same expression of Oi, O2 annihilation contribu- 
tions. O5, Oe have & {V — A) (gi {V + A) vertex so they have 
chiral enhancement compared to {V — A) ^ {V — A) vertex, 
and its' CKM factor are V*i,Vts, then its' contributions are 
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M. 



p(d) 



- -F^'k, 



sVeQq/sTi" 



4M| 



-Tk* / dX2 j 



b2db2 



X [{1 - X2Wk4^2) - {1 + X2Wk4^2)] 



t„ = max 



{Da,l/b2)) 



(60) 



(c) 




(d) 




FIG. 11: Annihilation diagrams caused by O3 ~ Oe operators. 
QCD penguin annihilations include Os in Wilson coefficient, 
so they are the same order of all contributions except for Oi , 
O2 annihilations. 



comparatively large and get main origins for isospin breaking 
effects. 



j^O(a, ^ j^r 



tk' I dxi / bidb 



^ ^' 4M| 

(i::=max(yl„,l/6i)) (55) 



M. 



S(b) 



p(b) 



= -M, 

= -F^'kt^-^^§^l'dx2jb2db2aUtl) 

X St(x2)e-[^^*(*'']0j.(a;2)iJjf<^'(62Ba) 



ft' = max(Ba, l/&2)j 



(61) 



M, 



S(d) 



p(d) 



-M, 



-F(«)e 



3\/6Qq/s7r /■! 



2A/| 



' da;2 / b2db2 aeitt) 

X S,(a.2)e-[«-*(*^'V|^.(x2)zjH«(&27?.) 
uf = max(Da, 1/62)) 



Al^{l + xi)Ml, Bl^{l-X2)Ml, 



Cl = xiM|, 



D^ = a;2A4 



(62) 



(63) 



^^^4 " ^ ^* AMI 



TK' 



f dx2 I 



62^62 



04 



(t^)5,(x2)e[-^-(*^)]^^ff(^)(62B.) 



X [(2-x2)<^K-(a;2)+a;20K-(a:2)] (56) 



M, 



p{b) 



TK* I dx2 I b2db 



a4{t)St{x2)ei-^'<'^<'>k^H^^\b2Ba) 

X [X2(I)K'(X2) + (2- X2)(t>K'ix2)] 

(tl = max(Ba,l/&2) 



(57) 



M. 



S(c) 



-Af 



^(c) 



-^4 



SV&QqfK'TV 



AMI 



M, 



S(d) 



a4(t:)5't(a;i)e[-^«(*°'Vi3(2^i,fei)^o(biC„) 
(i;^ = max(Ca, 1/61)) (58) 



a4(i^)5,(x2)e[-«^*(*»']i|//(«(fe2^a) 
X [(l + a;2)0K->(2:2)-(l-X2)0K-(a;2)] (59) 



V. LONG DISTANCE CONTRIBUTIONS TO 
THE PHOTON QUARK COUPLING 




h -^ 



FIG. 12: Vector-Meson-Dominance contributions mediated 

by V'.P.^- 

Here we want to discuss the long distance contributions. In 
order to examine the standard model or search for new physics 
indirectly by comparing the experimental data with the values 
predicted within the standard model, we have to take into ac- 
count these long distance effects: B — > K*{J/tp,p,ij) — > K*^ 
[B, 18] (Figinj. It should be noted that B -^ DDK* ^ K'j 
is small compared to the J/tp intermediate state contribution. 

These contributions are caused by Oi, O2 operators, and 
the effective Hamiltonian describing these processes is 

H^ff = ^ E V,tV:,{C^it)0[''\t) + C2{t)0[^\t)) + h.c. 

(64) 

If we use the vector-meson-dominance, the B — > K*^ decay 
amplitude can be expressed as inserting the complete set of 
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(A) 



(B) 




u u 




(C) 



(D) 





FIG. 13: (A), (B) are factorizable and (C), (D) are nonfac- 
torizable contributions to tiie hadronic matrix element for 

<K*ij\H,ff\B>. 



possible intermediate vector meson states like 



(7^*7l^e//|B)=^{7l^.^e^!V) 



.^^(yK*\H^jf\B), 

HV '"'V 



(65) 



where V = ■(/), p,cij. Now we concretely consider the 

B -^ K*^ — > A'*7. Four diagrams contribute to the hadronic 
matrix element of (-ft'*V'|^e//|B) (see Fig |13ll . and first of all, 
we consider the leading contributions: the factorizable ones, 
Figs. ElA) andHSlB). 



A. Factorizable contribution 

In this case, the B -^ ipK* decay amplitude can be factor- 
ized as 

{i^K- I H^ff I B) = ^K,v;:ai(t){^ i C7^(l - 7^)c | 0) 
x{A:* 1 S7''(l-7'''')6| B), (66) 

and the definition of the decay constant is 

{V'(g) I C7^c j 0) = im^g^{q^)e'^^{q), (67) 

then the decay amplitude can be written as 
G 



M{B -^ K*'4>{q)) = -j=VcbV*sai{t)im^g^{q )e*^f,{q) 



V2 

x{K' i S7''(l-7^)b| B), 



(68) 



where a\{t) — Ci(t) + C2(i)/3. The conversion part of the i/) 
meson into photon can be expressed as 



(7 I A^J^,„ I tjj} = --em^g^{q'), 



(69) 



then the total amplitude of B — * K*^ mediated by tp meson 
can be expressed as follows, 

M(B ^ if > ^ 7^*7) = ^i/,y>i(t) (^H£»|(2)!^ 4^ 

x{K* \s'y''{l-j'')b\B) (70) 



where the real photon momentum is q^ = 0. In principle, we 
need to include the width of the vector meson in the propa- 
gator and write 



(71) 



g2 — m'-^ + imV 

but we have (F^/m^) ~ O(10~^) and the effects of the width 
can be safely neglected. 

The amplitudes for B — > K*ij -^ K*^ can be computed in 
a similar manner. In this case, {ru/niu) = 1.0 x lO"'^ and 
we can also neglect the width effect in the meson propagator. 
Differences with B -^ K*%l) -^ K*y are the value of decay con- 
stant g^ (0) and the factor for the electromagnetic interaction. 

M{B-^K'lo-^K*j) = ^KbKr.ai(i) ('^^^ 

X (A'* I 57^(1 -7^)6 I B) (72) 

However in the B —* K* p —* A'*7 case, the p resonance peak 
is not so sharp, so the propagation of p meson generates the 
strong phase: (Fp/mp) ~ 0.19 and it introduces ~ 11° strong 
phase. 



M{B^ K*p^ K*'y) 



^K.K.ar(i) |^2(l-»r,/m,) 
x{A* I S7''(l-7^)6| B) 



(73) 



In order to estimate these long distance contributions, we 
have to know the decay constant gv- The decay constants are 
experimentally determined by the V -^ e^e~ data [1|. The 
amplitude for V — > e^e~ can be expressed as 



M(V ■ 



}e^mvgv{q^) 



(74) 



where Q expresses the electric charge like that Q = Qc when 
V — ip, Q — {Qu — Qd)/V2 inV = p case, and inV = uj case, 
Q = (Qu + Qd)/V2. Then the decay width for V ^ e+e" 
decay can be written like 



r(V-^e+e") = 



4TvQ^al^gy(q^ 



3mv 
and the values of gv are in Table Q 



(75) 



V 


r{V ^e+e-){GeV) 


mv{GeV) 


gl{GeV^) 


J/^(15) 


5.26 X 10"^ 


3.097 


0.1642 


^{2S) 


2.19 X lO"** 


3.686 


0.0814 


ip (3770) 


0.26 X 10-** 


3.770 


0.0099 


1^(4040) 


0.75 X 10-8 


4.040 


0.0306 


V'(4160) 


0.77 X 10"'' 


4.160 


0.0323 


i/;(4415) 


0.47 X lO-** 


4.415 


0.0209 


P 


7.02 X 10^8 


0.771 


0.0485 


LU 


0.60 X lO-** 


0.783 


0.0379 



TABLE I: The coefficients gv- 
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Furthermore, these decay constants are defined at the q^ = my energy scale. We need ones at q'^ = 0, so we have to 
extrapolate these decay constants frorn q^ = m^ to q^ = 0. We express gv{0) as ov-fO) = Kgv{q^) by using suppression factor 
K. In the ip cases, we take k ~ 0.4 |17lll8|l . and in the p,uj cases, we take k ~ 1.0 p^ 120(1 . Then the long distance contributions 
mediated by tp,p,ij are 



MiB^K^-y) = ^a^{t)e[v.,K:^-^^^^^p^ + V^,V:, 



g^imir 



9pi-ml) 



2\2 
p) 



)eUK*\sj^il-j')b\B), (76) 



6 2(l-ir/mp)_ 

and if we calculate the form factor of (A"* | S7''(l — 7'')b | B), the long distance contributions become as follows: 



M 



S(A) 



= -M 



P(A) 



S^%(0) '' 



Ml 






XrK* [(J)"k'{x2) + (i)K'{x2)] i 



-iu 



g.{miY 



gpimp) 



6 2{l -tV/mp) 



(i? = max{A7,B7,l/bi,l/b2)) 



(77) 



M 



S(B) 



^f(°) / dxidx2 /'&id6i62d&2)<^s(a;i,bi)Sf(x2)a,(4)e'-^«('''-^«*<*'"ai(4)H^''^(yl76i,C76i,C762) 
-'"s Jo J 

X ^{X2 +2)rK-'(l>K*{x2) +(^K. (3^2) -X2rK'cl>K*{x2)^ I £,c g '^ 



+ (u 



gu.(mlY ^ gpim^Y 



2(1 -iV/mp) 



(78) 



M 



P(B) 



8n^ 
'Ml 



f'°' / dxidx2 / bidbib2db2<l>B{xi,bi)St{x2)as{t7)e 



b\J-SB(t\)~SK' (4)1 „, (fb\ rr(b) 



X - X2rK' (l>K* {X2) + (f>K* (3^2) + {X2 + 2)rK* (l>K* i^l)] C 

(tr = max{A7, Ct, I/61, 1/62)) 

I 



aj{4)H!,''>{Arb,,C7bi,C7b2) 



2Kg^(m^) 



+ (u 



6 



+ 



Opijnlf 



2{1 ~iV/mp) 



(79) 



B. Nonfactorizable contribution 

Next we estimate the the effect of nonfactorizable con- 
tributions to the physical quantity like branching ratio, CP 
asymmetry, and isospin breaking effects. In order to do so in 
the case of _B — > K*il) -^ K*'y at first, we use the experimen- 
tal data on the branching ratio and different helicity ampli- 
tudes for B -^ J ji^) K* decay mode. The branching ratio is 
Br{B° -* J/ip K'°) = (1.31 ± 0.07) x lO'^Q, and the frac- 
tion of the transversely polarized decay width to the total 
decay width is about Tt/T =~ 0.4 [Ulillg^l, then the cor- 
responding transversally polarized branching ratio amounts 
to 



BriB -> J/i) K*)t ~ 5.0 x 10" 



(80) 



On the other hand, if we compute the branching ratio by using 
ea. l66L we have 



Br{B -^ J/ip K*)t ~ 2.3 x 10" 



(81) 



If we assume that the difference between the experimental 
value ea. l80ll and our prediction ea. ll81t is due to the nonfac- 
torizable amplitude, then 



nonfactorizable A{B — > J/xj) K*)t 
factorizable A{B -^ Jf-ip K')t 



:0.4 



(82) 



Note however, that B -^ K*'y is dominated by the short dis- 
tance amplitudes. The long distance correction from the fac- 
torizable diagram is about 4% of the total decay amplitude. 
So when we add the nonfactorizable amplitude, the long dis- 
tance correction increases to 6% in the total amplitude, and 
12% in the branching ratio. We have included these correc- 
tions in our numerical estimates given below. 

Furthermore, we estimate the effect of nonfactorizable con- 
tribution to the direct CP asymmetry. In general, a nonfac- 
torized amplitudes has a relative strong phase compared to 
the factorized amplitude. We already know that the nonfac- 
torizable diagram amounts to about 2% to the short distance 
amplitude, then we can numerically estimate the CP asym- 
metry uncertainty from the nonfactorizable diagram by in- 
troducing the strong phase as a free parameter. We conclude 
that only less than 10% uncertainty is generated by the long 
distance nonfactorizable amplitude, and as we will see later, 
this error is small compared to the total uncertainty in CP 
asymmetry from other origins. Finally, we mention that these 
long distance contributions do not generate the isospin break- 
ing effect, the nonfactorizable contribution can be neglected 
in computing the isospin breaking effects. 

In the case of B ^ K*{p,uj) -^ K*'y, we can expect that 
the factorized amplitudes are dominant to the total decay 
amp litude in B ^ K* (p, uj) by the analogy oi B ^> pp decay 
[24|. then we can neglect the nonfactorized contribution to 
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(a) 




(b) 




C. Another diagrams for long distance 
contributions to the photon quark coupling 



FIG. 14: Long distance effects mediated by p,uj which con- 
tribute only to the charged mode. 



the above physical quantities. 



Next we want to consider another contribution with dif- 
ferent topology which exist only in the charged decay mode 
like B -^ K* J ('Fia ll4ll . If we neglect the nonfactorizable 
contributions and annihilation contributions, there are two 
diagrams that contribute to the hadronic matrix elements 
{K* p{lo) I Heff I B ). We define the p or oj meson momen- 
tum P3 = Mb/V2 (1, 0, Ot) and the spectator quark momen- 
tum fraction as x^. 



47r2/K- 



M 



2 -^ ? 



Jo J 11 -iT/mp 3 J 



X a2inrp[<t>;{x3)+<t>;{xa)] H^^^A^hi^B-jW^Bjbs) (t" = max(A7, B7, 1/61, 1/63)) 



(83) 



M 



S(.) _ ^^^^^(0)^^ Cdx^dx, [b,db,hdh <t>Bix,M)\ , ^^^^'P + g^jg.(x2)a.(t'')e[-^^''^>-^>'<''')] 
'0 J 11 — iT/mp 3 J 



M 



p(b) 



Mi 



X a2(t') [(a:3 + 2)rp,?i^(x3)+0j(x3)-a:3rp0^(:r3)] i:f^"(A7fei, Cyfei, C7&3 



(84) 



47rVK' 
M| 



■F'°^^u I dxidX3 



bidhibsdbs (l)B{xi,bi)\- —7 1 ^ \St{x2)as{t )ey ^' ' "^ 'i 

L i — zi /nip 6 J 



X a2 



it") [-a;3rp<^;(a;3) + 0p(a;3) + (xa + 2)rprp{xs)\ H^''\Arbi,C7bi,C7bs) 
(.b -_(A7,C7, 1/61, 1/63)) 



t = max( 



(85) 



A7 = XxX3M%, 



B? ^xiMl, 



C? = XiMl 



(86) 



In the computation of the above formulas, we use the p and u) meson wave function extracted from light-cone QCD sum rule 
im . and the detailed expression is in Appendix|n| 



VI. NUMERICAL RESULTS 

We want to show the numerical analysis in this section. In 
the evaluation of the various form factors and amplitudes, we 
adopt Gf ~ 1.16639 x 10~*GeV~^, leading order strong cou- 
pling fla defined at the fiavor number n/ = 4, the decay con- 
stants /fl = 190MeV, /k* = 226MeV, and fj^. = 185MeV, 
the masses Mb = 5.28GeV, Mk' = 0.892GeV and 
rric = 1.2GeV, the meson lifetime rgo = 1.542 ps and 
Tg+ — 1.674 ps. Furthermore we used the leading order 
Wilson coefficients 8] and we take the K* , p, and uj meson 



wave functions up to twist-3. In order to make clear the 
theoretical error of the predicted physical quantities, we want 
to show how to estimate these errors. 

A. Error Estimation 

When we estimate the physical quantities like branching 
ratio, CP asymmetry, and isospin breaking effect, there are 
four major classes of error in pQCD computations: (1) the 
input parameter uncertainties; (2) higher order effects in per- 
turbation expansion; (3) the CKM parameter uncertainties; 
and (4) the hadronic uncertainties from the u quark loop. 



J 



1. First we want to estimate the class(l) error for various physical quantities. For class(l), we change the decay constants, 
the B meson wave function parameter ujb, and c parameter of the threshold function. We estimate the uncertainties from 
the decay constants to be 15% in the amplitude. If we change the ub in the range ujb = (0.40 ± 0.04) GeV, and c in the 
range c = 0.4 ± 0.1, these uncertainties change the B — > K* form factor by about 15% at the amplitude level. Thus we 
regard the total uncertainty for class(l) to be 20%. Here we discuss how this error affects the experimental observables 
such as the branching ratio, direct CP asymmetry, and isospin breaking. 
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Branching Ratio 

In order to see how much error is generated when we change some parameters in class(l), we introduce real parameter 

61 's as the fractional differences of the amplitudes from ones with a fixed hadronic parameter, where i and j express 

the flavor and electric charge. Note that the uncertainty in decay constants leads to an uncertainty in overall factor 

of the amplitude, i.e. they don't lead to an uncertainty in the phase of the amplitude. In the change wave function 

parameters on the other hand, the phase changes a little, but it's effect is very small and we can introduce S^'s as 

real parameters. 

The decay widths of the B and B meson decays can be expressed as 

r{B') = \v;,VtsAl{l + Si) + v:,V.sAi{l + 5i)+v:,VnsAi{l + 5i)f (87) 

r(B^) = IVtbVtUlil + S',) + VctV:,Aiil + Si) + VubV:sAiil + 5i)\^ 



and we can see that the uncertainty to the branching ratio from input parameters comes from the error of the Or-y 
amplitude, and it amounts to about 25-[ ~ 40%. 

• Direct CP Asymmetry 

From eg. 11871 and 18811 . the direct CP asymmetry can be expressed as follows, 

AcP = ,^^ ,^ ,,, f,,^, — \lm{V:tVtsVcbV:,)ImiAlA:'){l + Si) + /m(l^,;KsKi,K*s)/m(A^^:^)(l + Si)] (89) 

and the error for it is 



AAcP _ A'cp - AcP ^^ {Sj ^ Si)Im{Ai^/Ar) + jSj - Si)ImiAi*/Ai*) 
Acp Acp '^ {l + 5l)[Im{Ai''/Al*) + Im{Ai*/A{'')\ 



(90) 



We can see that the uncertainties can cancel. We have checked that numerically the class(l) error for the CP 
asymmetry amounts to few percent and is small compared to other errors (see below). 

• Isospin Breaking 

On the other hand, we want to show that the hadronic parameter uncertainties especially from up and c dependences 
of the isospin breaking effect can be large even though we take the ratio as the CP asymmetry. The decay width 
of the neutral and charged decay modes with the theoretical error can be written from ea. (l87ll and I88II as 

r" = \V,lVts^l{l + S'l) + v:,V.s^l{l + 5'-l) + v:,VusAl{l + 5l)\\ (91) 

r+ = \VtlVtsAt{l + St) + v:bV.sAt{l + St) + V:t,VusAt{l + St)\\ (92) 

and the isospin breaking effect is given by 

w _ \A'A\^ + 5'i)'-\At\\l + 5t)' ,.., 

"+ j^?|2(l + 50)2 + |A+P(l + 5+)2 ' ^''^ 

where we neglected all terms except for those proportional to \Al\^ because |AJp/|^^p ~ 0(10"''), and the CKM 
factor of the |yli|^ is suppressed as IKibTCs/VibytsI^ ~ 0(A*). Then the error can be expressed as 



A(Ao+) _ A^, - Ao, ^ 4|A»|V+P[^°^^+] 



i0+ (iA?P-|4+12)(|A?|2 + iA+P) 



(94) 



We can easily imagine that the decay constant uncertainties are canceled as the direct CP asymmetry. How- 
ever we observe that even though (5° — Sf is small, there exist | A? p — | Af \ ^ in the denominator and it is 
also small, then the error enhancement can occur. Variation of lob and c introduces S'^ — Sf ~ 0.5% while 
{\A°\^ - \At\'^)/{\At\^ + l^^l^) ~ 5%. This gives about 20% error for the isospin breaking. From the above argu- 
ment, we can see that the error from LOp and c uncertainties remain somewhat large. Thus we estimate the class(l) 
error for the isospin breaking effect to be about 20%. 

Next we want to discuss the class(2) error. For class(2), we expect an error coming from the fact that we used the 
leading order term in ois{t). There are also errors coming from neglecting higher order decay amplitudes. But we have 
not checked the effect of class(2) errors as it requires actual computation of higher order amplitudes. We guess that the 
error is approximately 15% in the amplitude. Then the theoretical errors from class(2) are 30% in the branching ratio, 
about a few % in the direct CP asymmetry, and 20% in the isospin breaking effect. 

About the class(3) error, we change the p, f) parameter in the range p — p(l — A^/2) = 0.20 ± 0.09 and 
fl — ■q{l ~ A'^/2) = 0.33 ± 0.05 01, and numerically estimate how the physical quantities are affected by the changing 
of parameters. The major contributions to the branching ratio and isospin breaking effects come from the terms which 
are proportional to Vj^Vts, so they are less sensitive to the error in p, fj. On the other hand, direct CP asymmetry 
depends on IraiV^tytsVcbV^g) and Imiy^iytsVubV^s) as in eg. 11891 . thus the error from the p, fj uncertainties amounts to 
about 15%. 
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4. The class(4) error comes from the u quark loop hadronic uncertainties. The terms which are proportional to V^i,Vus are 
not very important to the computation of the branching ratio and isospin breaking effect, so for these quantities we can 
neglect the class(4) uncertainties. 

However for CP asymmetry, c and u quark loops give comparable contributions as seen in eg. 18911 . thus the u quark 
loop contribution which is infected with nonperturbative correction, cannot be neglected. If we regard the u quark loop 
uncertainty as about 100% at the amplitude level for both real and imaginary part, the numerical error for the direct 
CP asymmetry amounts to about 75%. 

In summary, we regard the error of the branching ratio, direct CP asymmetry, and isospin breaking effects as 50% (class(l); 
40%, class(2); 30%), 75% (class(4); 75%), and 30% (class(l); 20%, class(2); 20%), respectively. 

I 



KftKi 



FIG. 15: CKM unitary triangle 



B. Numerical Results 

The numerical results for each decay amplitude Mi in the 
neutral decay (Tabl^HJ and charged decay (Tabl JTTTl in unit 
of 10"'^GeV"^ are as follows. 

The total decay amplitude can be expressed by using these 
components as 



A{B^ K*-/) 



+ 



Mt + Mc + Mu 

(e; ■ e*j,,)(Mf + M^ + M^ 

ie 



^,+-e;^e'i^, (Mf + Aff + M^ 



(95) 



Mt 

Me 

Mu 



Mr-y+Msg + Ma^e, 

Mlu+2u + M2+Mp+^, 



where all components include CKM factors. If we express 
K* and 7 helicities as Ai,A2, the combinations which can 
contribute to the decay amplitude are j4ai,A2 = ^+,+,^-,-, 
if we take into account the fact that B meson is spinless and 
a real photon has helicities ±1. Then the total decay width 
of _B ^ K*'y is given by 



r = 



Mf + M! + M^\^ 



8tvMb 
and the branching ratios for B 



IMf + M^ + M^f 



,(96) 



K*"f become as follows: 



Br{B° ^ K'^'y) = (5.8 ± 2.9) x 10~^ (97) 

Br{B^ -^ K*'^-y) = (6.0 ± 3.0) x 10"^ (98) 

Next we want to extract the direct CP asymmetry. We take 
into account up to 0{X'^) about the CKM matrix components. 



Vkm — 




f 1 - AV2 - 


-^ A AX\p-tr,)\ 


-A 


1 - AV2 - (1/8 + AV2) A* AX' 


\a\\i^p- 


-iri) -AX' + AX'' {1/2 -p-iri) l-A'X''/2J 




(99) 

1 



and the unitary triangle related to this decay mode should be 
crushed ('Fig ll5|l . 

If we express each amplitudes Mi as £,iAie'^' where 
?i = Vi'bVis/V*i,Vcs, in order to separate weak phase and 
strong phase Si, the decay amplitudes can be rewritten as 



A{B 
A{B 



K*'y) = V:tV,s [CtAte"* +^cA,e'''' 



'S,uAue" 



K*l) = V^bV:, [CtAte''' + CcAce''^ + C^^e'*" 



and the direct CP asymmetry can be expressed as 

r(g -^ K*j) - TjB -^ K*j) 



Acp 



r{B -^ K-'-i) + r(B ^ K*'y) 



Rn 
Rd' 



(100) 
(101) 



(102) 



Rt. 



Rd 



\^AtAc sm{5t - 5,)Im{VtbVt:v:tVcs) 

+A,Au sm{5c - 5n)lm{VcbV:,V:tVus] 

+AuAt sin((5„ - 5t)lm{VnbVZsVttVts) 

{A'tlVtbVtlf + Al\V,bV:f + Al\V^bV:,\ 
+AtA,cos{5t ~ S,)Re{VtbVtlV,lVcs) 
+AcAu cos(5c - Su)Re{VcbV*sVubVus 
+AnAt cos(<5„ - St)Re{VubV:,V,lVcs] 



then its' values are 



Acp{B° ^ K*°-y) = -(6.1 ± 4.6) X 10" 
Acp{B^ ^ K*^j) = -(5.7 ± 4.3) X 10" 



(103) 

/2 

(104) 



(105) 
(106) 



Finally, we want to estimate the isospin breaking effects 
as eq.(|5J. This effect is caused by Osg (Fig^J, c and u loop 
contributions (Fig|^, Oi ~ Oa annihilation fFigs llO| and llH . 
and the long distance contributions mediated p and ui in 
charged mode (Fig |14ll . About the bremsstrahlung photon 
contributions emitted through quark lines, whether the spec- 
tator quark is it or d affects the strength and the sign for 
the coupling of photon and quark line, so they generate the 
isospin breaking effects. The most important contributions to 
the isospin breaking effects come from QCD penguin O5, Oq 
annihilation. These effects are additive to the dominant con- 
tribution O7-, in both neutral and charged decays (see Tables 
HTIand llTB . However its' size are different: the neutral mode's 
is larger than charged mode's. Then the sign of total isospin 
breaking effects becomes plus and it's value is as follows. 



Ao+ = +(2.7 ±0.8) X 10" 



(107) 
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Mf/F^''> 


Mf/7^(«) 


VtWu 


MfJF^°'> 


m£,/f(«) 


ACe/i^*'" 


Aff^/Ffo) 


Af|;/F(«) 


A'Ce/^*"' 




-218.67-3.86i 


-2.19-0.55i 


-ll.56-5.70i 


218.67+3.86i 


2.27-H0.59i 


11.58-H5.63i 


V.lVcs 


MfjF'^°'> 


MijF^°'> 


M^/F^°^ 


A/f,/F(°) 


Mi'jF'-°'> 


a^^/f(°) 




-0.29-l.Oli 


6.42-12.63i 


-13.29 


-0.19+1. 27i 


-4.81-l-8.23i 


15.09 


v:,Vus 


MiL+2„/f(0) 


Mi/F^°^ 


A./;+./i^'°) 


aC+2„/f(») 


Mi/F^"'' 


Afp^+./F(«) 




-0.63-K0.22i 





-0.03-0.06i 


0.67-0.18i 





0.03-f0.07i 




TABLE II 


. 50 ^ ;^.0 


7 at p = 0.20, r; = 0.33, lob = 0.40GeV. 




Mf/F^"^ 


A/f/F(°) 


VtlVts 


M7^^/f(°) 


a4/f<"' 


ACe/i^*"' 


Aff^/F(°) 


A/i;/i^<°' 


mL,/f'''^ 




-218.67-3.86i 


-4.89-O.lOi 


-2.47+0.37i 


218.67+3.86i 


4.83-0.82i 


2.86-h0.14i 


v:,Vcs 


A.f£/i^(°) 


Mf,/F(«) 


mS/f^'^ 


A/fc/F<"' 


MijF^''^ 


A^;/F™ 




-0.66+2.15i 


6.42-12.631 


-13.29 


1.39-2.60i 


-4.8H-8.23i 


15.09 


Vy_l,VuS 


M,U,JF^°^ 


Af#/i^(°) 


A/p^+./F(»' 


aC+2„/f(°) 


A^r/F(°) 


a//j;./fw 




-0.75+0.51i 


0.35-H.Oli 


-0.04-f0.05i 


0.79-0.18i 


-0.75-fl.l6i 


0.05-0.05i 



TABLE III: B+ -^ K*+^ at p = 0.20, f) = 0.33, ub = 0.40GeV. 



VII. CONCLUSION 

In tiiis paper, we calculated the branching ratio, direct CP 
asymmetry, and isospin breaking effect within the standard 
model using the pQCD approach. It is useful to compare 
our results with those existing in the literature. The decay 
amplitude can be obtained from the transition form factor 

= -iTf {O)e,,o.ppefi. P^q" (108) 

where P = Pi -I- P2 , q = Pi — P2- Within the framework of 
pQCD, we obtain the value of the B —> K* transition form 
factor as Tf' (0) = 0.23 ± 0.06. The resuh can be com- 
pared with the ones extracted by another estimation. In 
the QCD Factorization, T/^* (0) = 0.27 ± 0.04 25] and an 
updated phenomenological estimate of this quantity with 
the light-cone distribution amplitudes for the K* meson 
is Ti^ (0) =0.27 ±0.02 ;26|. While the central value in 
the updated result is the same as before, the error is re- 
duced by a factor of 2. In the light-cone QCD sum 



rule Tf^ (0) 
Tr(0) 



0.38 ±0.06 |27|], the lattice QCD simulation 
+0.04 r^ „„j rr.K* /r,\ _ n oc;-l-0.05 p^ ^^^ ^^le 



0.321^:^* |2g and T{ 



■(0) = 0.251^2 11 



covariant light-front approach T/^ (0) = 0.24 [Sg. There 
are several estimates of the branching ratio by using the 
value of T^' (0) = 0.38 ± 0.06 extracted from light-cone QCD 
sum rule. Comparing the results with experiments, this 
value of the form factor over estimates the branching ra- 
tios IM El E3- Also, it should be noted that Tf (q^) 
and other related form factors have been computed in the 
frame work of pQCD |33 |. They obtained the central value 
as T^ (0) = 0.315. The difference between our results and 
their's is the K* meson wave function. We take the new K* 
wave function parameters computed in Ref . |l3| . 

Note that we have also included the long dis- 
tance contributions. If we neglect them, the branch- 
ing ratios become Br{B° -^ K*°"f) = (5.2 ± 2.6) x 10"^ and 



Br{B^ -* ^"^7) = (5.3 ± 2.7) x lO"'^ to be pared with re- 
suhs shown in Eq.linZJand||nHJ- The B ^ K*i> ^ K'j con- 
tribution to the total decay width amounts to about 12% and 
also it works additive to the branching ratios. We also em- 
phasize that we can calculate the annihilation contributions 
with the pQCD approach, and these contribute to the total 
decay width which amount to about 2 ~ 10%. 

This analysis predicts less than 1% direct CP asym- 
metry within the standard model. If we neglect 
the long distance contributions, the asymmetries be- 
come Acp{B° ~* K*°^) = -(6.7 ± 5.0) x 10"^ and 
Acp{B^ -^ K*^-)) = -(7.2 ± 5.4) x 10~^ and as to 
isospin breaking effect like Ao+ = ±(2.6 ±0.8) x 10"^. The 
long distance contributions do not seem to affect to these 
asymmetries. 

The branching ratio of the neutral decay is similar to that 
of the charged decay, in spite of the difference of the lifetime 
between them. This effect is mainly caused by the 4-quark 
penguin operators O5, Os- If we neglect these contributions, 
the isospin breaking is Ao+ = —(1.2 ± 0.4) x 10~^, so we can 
see that they generate about 4% isospin breaking effect. This 
result is similar to the conclusion of Ref. |35(| . 

B — > A'*7 decay, as we have first mentioned, is an attractive 
decay mode to test the standard model and search for new 
physics. In order to look for the new physics, we have to 
reduce the experimental errors. The error to the direct CP 
asymmetry must get smaller than 1%. That is to say, we need 
at least 20 times more data. This is not possible without the 
super B factory. 
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(a) 



(d) 




FIG. 16: 0{as) corrections to the hard scattering H . 
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FIG. 17: An electron which is scattered by the electromag- 
netic interaction (a) is observed with many soft photons. Sim- 
ilarly, a quark which is scattered by the strong interaction (b) 
is not observed as a single gluon: accompanied by many soft 
gluons, and they form hadron jets. 



APPENDIX A: BRIEF REVIEW OF PQCD 
1. Divergences in perturbative diagrams 

Here we want to review the kr factorization [3bI . At higher 
order, infinitely many gluon exchanges must be considered. In 
order to understand the factorization procedure, we refer to 
the diagrams of Fig ll6l 

They describe the 0{as) radiative corrections to the hard 
scattering process H. In general, individual higher order 
diagrams have two types of infrared divergences: soft and 
coUinear. Soft divergence comes from the region of a loop 
momentum where all it's momentum components in the light- 
cone coordinate vanish: 



r = (;+,/-, /t) = (A, A,A). 



(Al) 



CoUinear divergence originates from the gluon momentum re- 
gion which is parallel to the massless quark momentum. 



r = {i+,r,iT) ~ {Mb,a^/Mb,a). 



(A2) 



In both cases, the loop integration correspond to 
J d I l/l ~ log A, so logarithmic divergences are generated. 
It has been shown order by order in perturbation theory that 
these divergences can be separated from hard kernel and ab- 
sorbed into meson wave functions using eikonal approxima- 
tion [33. 

Furthermore, there are also double logarithm divergences 
in Fig |16f a) and llGr b*) when soft and coUinear momentum 
overlap. These large double logarithm can be summed by 
using renormalization group equation. This factor is called 
the Sudakov factor and also factorized into the definition of 
meson wave function pa . l39ll4Cl . The explicit expression for 
Sudakov factor is given by [331 (see Appendix [51. 

There are also ultraviolet divergences, and also another 
type of double logarithm which comes from the loop correc- 
tion for the weak decay vertex correction. These double log- 
arithm can also be factored out from hard part and grouped 
into the quark jet function. These double logarithms also 



should be resumed as the threshold factor |411l42l|. This fac- 
tor decreases faster than any other power of a; as a; —> 0, so 
it removes the endpoint singularity. Thus we can factor out 
the Sudakov factor, the threshold factor, and the ultravio- 
let divergences from hard part and grouped into meson wave 
function (Appendix |bJ . Then the redefinition of wave func- 
tions including these loop corrections get factorization energy 
scale dependence t. 

Thus the amplitude can be factorized into a perturbative 
part including a hard gluon exchange, and a nonperturba- 
tive part characterized by the meson distribution amplitudes. 
Then the total decay amplitude can be expressed as the con- 
volution: 

f dxidx2 d^hid^b2C{t)®^K*{x2,b2,t) 

Jo 

'SiH{xi,X2,bi,b2,t) 'Si $B{xi,bi,t), (A3) 

here ^K'{x2,b2,t), "I>s(a;i, 6i,i) are meson distribution am- 
plitudes that contain the soft divergences which come from 
quantum correction and H{xi,X2,b\,b2,i) is the hard kernel 
including finite piece of quantum correction, where 61, 62 are 
the conjugate variables to transverse momentum, and xi, X2 
are the momentum fractions of spectator quarks. 



2. Physical interpretation of Sudakov factor 

In order to understand the Sudakov factor physically, first 
we consider QED. When a charged particle is accelerated, in- 
finitely many photons must be emitted by the bremsstrahlung 
('Fia llTf a')'). A similar phenomenon occurs when a quark is ac- 
celerated: infinitely many gluons must be emitted. According 
to the feature of strong interaction, gluons cannot exist freely, 
so hadronic jet is produced. Then we observe many hadrons 
in the end if gluonic bremsstrahlung occurs. Thus the ampli- 
tude for an exclusive decay B — > K*"f is proportional to the 
probability that no bremsstrahlung gluon is emitted. This is 
the Sudakov factor and it is depicted in Fig |19l As seen in 
Fig |19l the Sudakov factor is large for small 6 and Q. Large 6 
implies that the quark and antiquark pair is separated, which 
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FIG. 18; b is the transverse interval between the quark and 
antiquark pair in the B meson. 




FIG. 19: The dependence of the Sudakov factor exp[—s{Q,b)] 
on Q and b where Q is the b quark momentum, and b is the 
interval between quarks which form hadrons. It is clear that 
the large b and Q region is suppressed. 



in turn implies less color shielding (see Fig |181 . Similar ab- 
sence of shielding occurs when b quark carries most of the 
momentum while the momentum fraction of spectator quark 
X in the B meson is small. 

Then the Sudakov factor suppresses the long distance con- 
tributions for the decay process a nd g ives the effective cutoff 
about the transverse direction |40ll43|. In short, the Sudakov 
factor corresponds to the probability for emitting no pho- 
tons. According to this factor, the property of short distance 
is guaranteed. 



APPENDIX B: SOME FUNCTIONS 



The expressions for some functions are presented in this 
appendix. In our numerical calculation, we use the leading 
order Oa formula. 



.(m) = 



27r 



/3o ln(/i/A„^. ) 



f3, . ^1^ (Bl) 



The explicit expression for Sudakov factor s{t,b) is given by 



s{t,b)=r ^ 

Ji/b M 



ln(^^jyl(a.(/.))+B(a,(M)) 



(B2) 



B 



Cf 


Qs 


+ 


/Qs^ 


2 


n 




\ TT J 




X 


'67 




2 


10 


9 




3 


27' 


2a,, 


J 


g27£; 


-1 



2. , /4 



3 TT 



(B3) 
(B4) 



where 7_b = 0.5722 is Euler constant and Cf ~ 4/3 is color 
factor. The meson wave function including summation factor 
has energy dependence 

<^s(3;i,6i,t) = </>s(a;i,fei)exp[-5's(i)], (B5) 

(j>K*(x2,t) = (l>K*{x2)exp[-SK*{t)], (B6) 

and the total functions including the Sudakov factor and the 
ultraviolet divergences are 

Seit) = sixiPr,bi) + 2 f ^7(q,(^)), (B7) 

Jl/bi M 
SK'{t) = s{x2Pi,b2)+s{{l-X2)P2,b2) 

+2 r ^7(«.(m)). (B8) 

Jl/b2 /^ 

Threshold factor is expressed as below [4ll lii]]. and we take 
the value c = 0.4. 

2i+2T(3/2 + c) 



St{x) 



0Fr(l- 



[x{l^x)r 



(B9) 



APPENDIX C: WAVE FUNCTIONS 

For the B meson wave function, we adopt the model 
1 



/2K 



i>B{xi,bi) — / dki d ki±e 



+ MB)Y(t'B{ki), 



(CI) 



4>B{ki) 



= A'^sa;i(l — x\) exp 



1 f xiMB 
'2\ LUB 



7,2 2 



(C2) 



with th shape parameter ujb ~ (0.40 ± 0.04)GeV. The nor- 
malization constant Nb is fixed by the decay constant /b 



dx<j)B{x, b = 0) 



fB 



2V2iVe 



(C3) 



where Nc is the color number. 

We use the vector meson wave functions determined by 
the light-cone QCD sum rule [VR H^. We choose the vector 
meson momentum P moving in the "-" direction along the z 
axis with P = My, and the polarization vectors el, et are 
defined as 



£i = (0,1,0), eT= (0,0,^(±l,-i)j, (C4) 

and er satisfies the gauge invariant condition P ■ er ~ 0. The 
nonlocal matrix elements sandwiched between the vacuum 
and the K* meson state can be expressed as follows. 



{K*-{P)\s{z)Iu{0)\0} = ^f^, 



x^^^^Mk. [ dxe 
P-2 Jo 



ixP-z O r (s) / \ 



(C5) 
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X / dxe^ '^(/)||(2;) + etm / dxe^ '^g\^ (x) 
Jo Jo 



(C6) 



{K*-{P) I s{z)to1M0) I 0) = -^/a" 



My. 



e^„pcy erPz I dxe -^9_l\^) 



(C7) 



{K'-{P)\ s(2)a^,it(0) |0) = -i 









(C8) 



where we neglect the terms proportional to r^* (twist-4) and 
the terms {niu + ms)/MK* ■ Then the K* meson distribution 
amplitudes up to twist-3 are 



^K'{P:eL) = ^ 



dxe"''-'{MK'[iL]^K'{.x) 



V^KJo 



(C9) 



$|^.(P,eT) = ^ 



2NcJo 
Mk' 
P-z 



dxe"' '''{Mk* [/t]0k* i^) 



+ [/t f]<^^.W +^ie^.p.[7V]eTP"^VK.W} 



(CIO) 



<j>K'{x) 



<t>K' (x) 



<t>K'{x) 



I 



K* 



2V2iVc 



^ ' 2y/2K " 



Ik' d ,(s) ,T / N Jk 

:«ll , <Pif. (a;) = 



AJWcdx 



2V27Vc 



= 31, <t>K'{x)^ 



2V2]Vc 



Ik' d (a) 



9>j2Kdx^ 



(Cll) 



where we use eoi23 = 1 and set the normalization condition 



{V) U) ,(t) ,(3)i 



about (l>i = {(l>n,(l}±,gY',g]'',h\ ',h\r} 



/ dx(j)i{x) 
Jo 



= 1. 



(C12) 



F 


A-* 


P 


MMeY] 


226 ± 28 


198 ±7 


fv[MeY] 


185 ± 10 


160 ± 10 


af [MeV] 


-0.4 ±0.2 





al|[MeV] 


0.09 ± 0.05 


0.18 ±0.10 


ailMeV] 


-0.34 ±0.18 





ai[MeV] 


0.13 ±0.09 


0.2 ±0.1 


s+ 


0.24 





S- 


-0.24 





h 


0.16 





6- 


-0.16 





S,3 


0.032 


0.032 


S,3 


0.013 


0.013 


<,3 


0.024 


0.024 


'^I'^O 


-2.1 


-2.1 



TABLE IV: Some parameter quantities. 



/i[''(a;) = 6a;(l-3;) 



1 + aiXi 



1 ± ,35 T 1 , 2 ,N 
ja2 + -g-C3 ) (5Si - 1] 



+35+ [3a::(l - x) +xlna::± (1 - x) ln(l - x)] 
+35- [a; In a:- (1 - i) ln(l - i)] 



h\:\a 



q2,3x /n2 i\|3x 2/r- 2 o\ 

Sxi + -ai Xi(3xi - 1) + -02 Xi (5xi - 3) 



sf W = 



+ ^Cr(3 - 30x? + 354) + h+ 1 + ^an , ^ 
4 A V i — a:: 

3 

+ -5-Xi [2 + Inx + ln(l - 2:)] 



6^(1 - :r) [1 + af e + {^4 + ^C, (l - ^^to 

+ ^C3}i5x'^i - 1)1 +6<5+[3a;(l-2:) + a;ln3; 

+ (l-x)ln(l-a;)] +65- [xlnx- (1 - x) ln(l - x)] 



g^l\x) = ^(l + 4)+af^4+0a^+5C3^)(34-l) 

± (y^4 + ^^3^ - f^^^-^to ] (3 - 304 ± 354 

3 - 

+ -5+ [2 + lux + ln(l - x)] 

3 - 

+ -5- [2xi + ln(l - x) - In i] 



I (a;) — 6a;(l — x) 



_(x) — 6x(l — : 



l + 3afx, + ^a!](5x?-l) 



l + 3aiXi + -a2 (5x- -1; 



Here Xi = 1 — 2x, and the expressions about p and a; me- 
son wave functions are the same as above with the values of 
parameters as follows evaluated at fi =lGeV (Tab liVt . Since 
p/lu states are {\uu) + |dd))/\/2, the qq distribution where 
q = u or d can be taken to the same for \uu) and \dd) using 
isospin symmetry. 



20 



S. Eidelman ei al. [Particle Data Group], Phys. Lett. B 
592 (2004) 1. 



[BELLE 
[BABAR 



Collaboration] 



Collaboration] 



[1] 

[2] K. Abe et al. 

arXiv:liep-ex/0409049 
[3] B. Aubert et al. 

,arXiv:liep-ex/0 408072 
[4] M. Nakao et al. [BELLE Colaboration Collaboration], 

Phys. Rev. D 69, 112001 (2004) arXiv:liep-ex/0402042 . 
[5] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

70, 112006 (2004) arXiv:hep-ex/0407003 . 
[6] J. Alexander et al. [Heavy Flavor Averaging Group 

(HFAG) Collaboration], arXiv:hep-ex/0412073 
[7] A. J. Buras, A. Czarnecki, M. Misiak and J. Urban, Nucl. 

Phys. B 631, 219 (2002) arXiv:hep-ph/0203135 . 
[8] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. 

Mod. Phys. 68, 1125 (1996) arXiv:hep-ph/9512380 . 
[9] M. Bauer and M. Wirbel, Z. Phys. C 42, 671 (1989). 
[10] Y. Y. Keum, H. n. Li and A. L Sanda, Phys. Lett. B 504, 

6 (2001) arXiv:hep-ph/0004004 . 
[11] P. Ball, V. M. Braun, Y. Koike and K. Tanaka, Nucl. 

Phys. B 529, 323 (1998) arXiv:hep-ph/9802299 . 
[12] P. BaU and M. Boglione, Phys. Rev. D 68, 094006 (2003) 

arXiv:hep-ph/0307337 . 
[13] M. Bander, D. Silverman and A. Soni, Phys. Rev. Lett. 

43, 242 (1979). 
[14] J. Liu and Y. P. Yao, Phys. Rev. D 42, 1485 (1990). 
[15] H. Simma and D. Wyler, Nucl. Phys. B 344, 283 (1990). 
[16] H. n. Li and G. L. Lin, Phys. Rev. D 60, 054001 (1999) 

arXiv:hep-ph/9812508 . 
[17] E. Golowich and S. Pakvasa, Phys. Rev. D 51, 1215 

(1995) arXiv:hep-ph/9408370 . 
[18] N. G. Deshpande, X. G. He and J. Trampetic, Phys. Lett. 

B 367, 362 (1996). 
[19] R. L. Anderson et al, Phys. Rev. Lett. 38, 263 (1977). 
[20] E. Paul, BONN-HE-81-26 Invited talk given at 10th Int. 

Symp. on Lepton and Photon Interactions at High En- 
ergy, Bonn, West Germany, Aug 24-29, 1981 
[21] K. Abe et al. [Belle Collaboration], Phys. Lett. B 538, 

11 (2002) arXiv:hep-ex/0205021 . 
[22] B. Aubert et al. [BABAR Collaborati on], Ph ys. Rev. 

Lett. 87, 241801 (2001) arXiv:hep-ex/0107049 . 
[23] M.Verderi, Workshop on the Unitarity Triangle in 2005, 

WG5 Session2, San Diego, 2005 (to be publushed) 
[24] A. L. Kagan, Phys. Lett. B 601, 151 (2004) 



[25; 
[26; 

[27 

[28; 

[29 

[30 

[31 

[32' 
[33 

[34; 

[35; 

[36 

[37 

[38; 

[39 
[40 
[41 

[42; 

[43; 

[44; 



'arXiv:hep-ph/0405134'. 

M. Beneke, T. Feldmann and D. Seide l, Nucl. Phys. B 

612, 25 (2001) arXiv:hep-ph/0106067". 

A. Ali, E. Lunghi and A. Y. Parkhomenko, Phys. Lett. 

B 595, 323 (2004) arXiv;hep-ph/0405075 . 

P. BaU and V. M. Braun, Phys. Rev. D 58, 094016 (1998) 

arXiv:hep-ph/9805422 . 

L. Del Debbio, J. M. Flynn, L. Lellouch and J. Nieves 
[UKQCD Collaboration], Phys. Lett. B 416, 392 (1998) 

arXiv;hep-lat/9708008 . 

D.Becirevic, talk presented at the Flavor Physics & CP 
Violation conference, Ecole Polytechnique, Paris France, 
June 2003 

H. Y. Cheng and C. K. Chua Phys. Rev. D 69, 094007 
(2004) arXiv:hep-ph/0401141'. 

A. Ah and A. Y. Parkhomenko, Eur. Phys. J. C 23, 89 
(2002) arXiv:hep-ph/0105302 . 
S. W. Bosch, arXiv:h ep-ph/0208203l 
S. W. Bosc h and G. Buchalla, Nucl. Phys. B 621, 459 
(2002) arXiv:hep-ph/0106081 . 

C. H. Chen and C. Q. Geng, Nucl. Phys. B 636, 338 
(2002) arXiv:hep-ph/0203003 . 
A. L. Kagan and M. Neubert, Phys. Lett. B 539, 227 

(2002) arXiv:hep-ph/0110078 . 

M. Nagashima and H. n. Li, Phys. Rev. D 67, 034001 

(2003) arXiv:hep-ph/0210173 . 

H. n. Li and H. L. Yu, Phys. Rev. D 53, 2480 (1996) 

arXiv:hep-ph/9411308 . 

J. C. Collins and D. E. Soper, Nucl. Phys. B 193, 381 
(1981). 

J. Botts and G. Sterman, Nucl. Phys. B 325, 62 (1989). 
H. n. Li and G. Sterman, Nucl. Phys. B 381, 129 (1992). 
H. n. Li, Phys. R ev. D 66, 094010 (2002) 

arXiv:hep-ph/0102013'. 

T. Kurimoto, H. n. Li and A. I. Sanda, Phys. Rev. D 65, 
014007 (2002) arXiv:hep-ph/0105003 . 
J. C. Collins and G. Sterman, Nucl. Phys. B 185 (1981) 
172. 

Y. Y. Keum, H. N. Li and A. L Sanda, Phys. Rev. D 63 
(2001) 054008 arXiv:hep-ph/0004173 . 



